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Abstract. Since the baryon-to-photon ratio 7710 is in some doubt at present, we ignore 
, the constraints on r)io from big bang nucleosynthesis (BBN) and fit the three key cosmo- 

. logical parameters (ft, nM,»?io) to four other observational constraints: Hubble parameter 

{ho), age of the universe (to), cluster gas (baryon) fraction (/o = fah^^^), and effective 
\l ' shape parameter (Fo). We consider open and flat CDM models and flat ACDM models, 

testing goodness of fit and drawing confidence regions by the Ax^ method. CDM models 
with Q,M = 1 (SCDM models) are accepted only because we allow a large error on ho, 
permitting h < 0.5. Open CDM models are accepted only for SIm > 0.4. ACDM models 
give similar results. In all of these models, large 7710 (> 6) is favored strongly over small 
riioi<y 2), supporting reports of low deuterium abundances on some QSO lines of sight, 
, and suggesting that observational determinations of primordial ^Ife may be contaminated 

' by systematic errors. Only if we drop the crucial Fo constraint are much lower values of 

' JIm and 7710 permitted. 

00 ! 

^ Key words: Baryon-to-photon ratio. Universal matter density, Hubble parameter 

Q ■ 1. Introduction 

■ In the context of the hot big bang cosmology, if the number of hght-neutrino 

species has its standard value Ni, = 3, the predicted primordial abundances 
of four light nuclides (D, ^He, ^He, and ^Li) depend on only one free param- 
^ . eter, r/10, the universal ratio (at present) of nucleons (baryons) to photons 

^ I (in units 10^^^). In principle, rjiQ is overdetermined by the observed or 

inferred primordial abundances of the four light nuclides. Indeed, Steigman, 
Schramm, and Gunn (1977) have exploited this fact to use big bang nucleo- 
synthesis (BBN) to constrain N^. The status quo ante is that observations, 
principally of D and ^He, have rendered r/10 one of the best known of the 
key cosmological parameters: rjiQ = 3.4 ±0.3 (Walker et al, 1991; the error 
bars being roughly "Icr"). 

At present, when the microwave background temperature T = 2.728 K 
(Fixsen et al, 1996), the universal baryonic mass-density parameter I^b ( = 
SttCpb/'^Hq ) is related to r/10 by 

nB = 3.667 X 10"^ r/10 = 0.0125 ± 0.0011. (1) 
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However, recently there have emerged reasons to suspect that rjiQ may not 
be so well determined, and even that the standard theory of BBN may 
not provide a very good fit to the current data (Hata et al., 1995). There 
are several options available for resolving this apparent conflict between 
theory and observation. Although some change in standard physics could 
offer resolution (e.g., a reduction in the effective value of N^, during BBN 
below its standard value 3; cf. Hata et al., 1995), Hata et al. (1995) note 
that large systematic errors may compromise the abundance data (cf. Copi, 
Schramm, and Turner, 1995). 

This controversy has been sharpened by new observations giving the deu- 
terium abundances on various lines of sight to high-redshift QSOs. In princi- 
ple, these data should yield the primordial D abundance, but current results 
span an order of magnitude. If the low value (D/H by number ~ 2 x 10~^; 
Tytler, Fan. and Buries, 1996; Buries and Tytler, 1996) is correct, then 
?7io ~ 7 in the standard model, but then it seems impossible to reconcile 
the inferred abundance of ^He [Olive and Steigman, 1995 (OS)] with (stan- 
dard) BBN for this value of r/io unless there are large systematic errors in 
the ^He data. If, instead, the high figures (D/H « 2 x 10~^; Carswell et 
al., 1994; Songaila et al., 1994; Rugers and Hogan, 1996) are correct, then 
D and ^He are consistent with r/io ~ 2, but modellers of Galactic chemi- 
cal evolution have a major puzzle: How has the Galaxy reduced D from its 
high primordial value to its present (local) low value without producing too 
much ^He (Steigman and Tosi, 1995), without using up too much interstel- 
lar gas (Edmunds, 1994; Prantzos, 1996), and without overproducing heavy 
elements (cf. Tosi, 1996, and references therein)? It appears that r/io, though 
known to order of magnitude, may now be among the less well-known cos- 
mological parameters. Despite this, large modern simulations which explore 
other cosmological parameters are often limited to a single value of rjiQ = 3.4 
(e.g., Borgani et al., 1997). 

Given this unsettled situation Steigman, Hata, and Felten (1997; here- 
after SHE) have proposed that it may be constructive to abandon nucle- 
osynthetic constraints on 7710 entirely and to put r/10 onto the same footing 
as the other cosmological free parameters, applying joint constraints on all 
these parameters based on other (non-BBN) astronomical observations and 
on theory and simulation. Armed with 7710 determined in this manner we 
may then "predict" the primordial abundances of the light nuclides and 
compare with the data to test the consistency of standard BBN. In this con- 
tribution to the proceedings of the ISSI Workshop on Primordial Nuclei and 
Their Galactic Evolution (Bern, Switzerland, 6-10 May 1997) we present 
a brief description of our approach along with a summary of our results 
for a "standard" (fiducial) choice of the observational constraints. For fur- 
ther details (especially of the many variations on the standard case to be 
described herein) and references the reader is encouraged to consult SHE. 
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Our approach is to let the three key cosmological parameters (/i, Om, "^/lo) 
range freely, fit the constraints (observables other than nucleosynthetic) to 
be described below, test goodness of fit by x^j and draw formal confidence 
regions for the parameters by the usual Ax^ method. Most of the SHF 
results are not surprising, and related work has been done before (White et 
al, 1996; Lineweaver et al, 1997; White and Silk, 1996; Bludman, 1997), but 
not with these three free variables and the full formalism. Further, some 
recent cosmological observations and simulations, particularly those related 
to the "shape parameter" T and the cluster baryon fraction (CBF), seem 
to pose a challenge to popular models, and there is some doubt whether 
any simple model presently fits all data well. Our approach, which begins 
by discarding nucleosynthetic constraints, provides a new way of looking 
at these problems. For example, the CBF and F constraints have not been 
applied jointly in earlier work which often also adopts a precise value for r/io 
(or rie). 

SHF find that, given their conservative (generous) choice of error bar on 
/i, the SCDM model is disfavored but by no means excluded. But even with 
this generous error bar, large values (^ 6) of ??io (J^b h"^ ^ 0.022) are favored 
strongly over low values 2; JIb ^ 0.007). This suggests that the low D 
abundances measured by Tytler et al. (1996) and by Buries and Tytler (1996) 
may be correct, and that the observed (extrapolated) primordial helium-4 
mass fraction [Yp k, 0.23; cf. OS and Olive, Skillman, and Steigman, 1997 
(OSS)], thought to be well determined, may be systematically too low for 
unknown reasons. 



3. CDM Models: Parameters and Observables 
3.1. Parameters 

The CDM models we consider are defined by three free parameters: Hubble 
parameter h; mass-density parameter f^M = SttCpm/^Hq; and baryon-to- 
photon ratio t/iq, related to Qb by equation (||). Here JIm by definition 
includes all "dynamical mass": mass which behaves dynamically like ordi- 
nary matter in the universal expansion; f^M is not limited to clustered mass 
only. Other free parameters having to do with structure formation, such as 
the tilt parameter n, could be added (White et al, 1996; Kolatt and Dekel, 
1997; White and Silk, 1996), but generally we have tried to avoid introducing 
many free parameters. 
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3.2. Observables 

We consider four observables (constraints) with measured values and with 
errors which are assumed to be Gaussian: (1) Hubble parameter ho] (2) age 
of the universe to; (3) gas-mass fraction /o = fch^^^ in rich clusters; and (4) 
"shape parameter" Fq from structure studies. In SHF we considered a fifth 
constraint: the dynamical mass-density parameter JIq a-s inferred from cluster 
measurements or from large-scale flows. Here, we ignore this constraint but 
we shall comment on its relation to our "standard" results. 

3.2.1. Observed Hubble Parameter ho 

For the Hubble parameter the observable is simply fit with the parameter 
h. Measurements of h still show scatter which is large compared with their 
formal error estimates (Bureau, Mould, and Stavcley-Smith, 1996; Tonry et 
al, 1997; Kundic et al, 1997; Tammann and Federspiel, 1997). This indicates 
systematic errors. To be conservative (permissive), we adopt ho = 0.70ib0.15. 
Perhaps a smaller error could be justified; below we will comment on the 
consequences of shrinking the error bar. 

3.2.2. Observed Age of the Universe to 

The age for these A = models is a function of h and Qm given by: t = 
9.78/1-1 x/(J)m;A=0) Qyj, [Weinberg, 1972, equations (15.3.11) & (15.3.20)]. 
We take the observed age of the oldest globular clusters as too = 14 it 2 Gyr 
(Bolte and Hogan, 1995; Jimenez, 1997; D'Antona, Caloi, and Mazzitelli, 
1997; Chaboyer et al, 1997; Cowan et al, 1997; cf. Nittler and Cowsik, 
1997). The universe is older than the oldest globular clusters by an unknown 
amount At. Although most theorists believe that At must be quite small 
(1 or 2 Gyr at most), we are unaware of any conclusive argument which 
guarantees this. To keep things simple, SHF introduced asymmetric error 
bars: to = lA'^l Gyr, allowing enough extra parameter space at large ages 
to accommodate a conservative range of At; extremely large ages will be 
eliminated by the ho constraint in any case. 

3.2.3. Observed Gas Mass Fraction in Clusters fo 

As is suggested by simulations, we assume that rich clusters provide a fairly 
unbiased sample of the universal ratio of baryonic to dark matter. 
Thus, we use the cluster (hot) gas fraction, /q, not as a constraint on Om, 
but as a constraint on the universal baryon fraction, the ratio Qb/^m- We 
emphasize that the following argument assumes that rich clusters provide a 
fair sample of the universal baryon fraction but does not assume that most 
of the mass in the universe, or any specific fraction of it, is in rich clusters. 

The measurement of /q poses some problems; for discussion and refer- 
ences, see SHF. We have followed the approach of Evrard, Metzler, and 
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Navarro (1996), who used gas-dynamical simulations to model the observa- 
tions. They find that the largest contribution to the error in /q arises from 
the measurement of the cluster's total mass, and they suggest that this error 
can be reduced by using an improved estimator and by restricting the mea- 
surement to regions of fairly high overdensity. Evrard (1997) applies these 
methods to data for real clusters and finds /g h^^"^ = 0.060 ± 0.003. To be 
conservative, we double his error bars and adopt for our constraint 

fo = /g h^^^ = 0.060 ± 0.006. (2) 

We are interested in the baryonic mass fraction, ^l^/^M, but even in 
rich clusters not all baryons are in the form of gas, and selection factors 
may operate in bringing baryons and dark matter into clusters. White et 
al. (1993) introduced a "baryon enhancement factor" T to describe these 
effects. T may be defined by 

fGo = rnG/^M, (3) 

where is the initial contribution of gas to 0,m (note that < JIb) and 
/go is the gas mass fraction in the cluster immediately after formation. T is 
really the gas enhancement factor, because the simulations do not distinguish 
between baryonic condensed objects if any (galaxies, stars, machos) and non- 
baryonic dark-matter particles. All of these are lumped together in the term 
(Om — ^g) and interact only by gravitation. 

If all the baryons start out as gas (JIg = ^b), and if gas turns into 
condensed objects only after cluster formation, then equation (^) may be 
rewritten: 

/G + /GAL = T17B/f^M, (4) 

where /g is the present cluster gas-mass fraction and /gal the present 
cluster mass fraction in baryonic condensed objects of all kinds (galaxies, 
stars, machos). White et al. (1993) took some pains to estimate the ratio 
/g//gal within the Abell radius of the Coma cluster, counting only galaxies 
(no stars or machos) in /gal- They obtained 

/g//gal = 5.5/i-3/2_ 

This is large, so unless systematic errors in this estimate are very large, 
the baryonic content of this cluster (at least) is dominated by the hot gas. 
Carrying /gal along as an indication of the size of the mean correction for 
all clusters, and solving equations (§) and (^) for fch^^'^, we find 

fch^/^ = [i + {h'/'/5.5)]-\TnB/nM)h'/', (6) 

where Qb is given from 7710 and h by equation (|^). 
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This is the appropriate theoretical function of the free parameters to 
fit to the observations. We set T = 0.9 in our "standard" case. Although 
this is representative of results from simulations, Cen (1997) finds that the 
determination of /g from X-ray observations may be biased toward high /g 
by large-scale projection effects; i.e., the calculated /g exceeds the true /g 
present in a cluster by a bias factor which can be as large as 1.4. Although 
Evrard et al. (1996) and Evrard (1997) have not observed such a bias in 
their simulations, SHF explored its effect on our analysis by using for T, 
instead of 0.9, an "effective value" T f« 0.9 x 1.4 ^ 1.3. 

3.2.4. Shape Parameter To from Large-Scale Structure 
The last observable we use is the "shape parameter" F, which describes the 
transfer function relating the initial perturbation spectrum Pi{k) oc A;" to 
the present spectrum P{k) of large-scale power fiuctuations, as observed, 
e.g., in the galaxy correlation function. When the spectral index n of Pi{k) 
has been chosen, T is determined by fitting the observed P{k). 

Results of observations may be cast in terms of an "effective shape param- 
eter" r (White et al., 1996) which we take as our observable. Studies show 
that for the usual range of CDM models, with or without A, the expression 
for r is 



(Peacock and Dodds, 1994; Sugiyama, 1995; Liddle et al., 1996a,b; White 
et al, 1996; Liddle and Viana, 1996; Peacock, 1997). For n 1, if 17b 
and r^B/l^M are small, we have F ~ JIm^- The Harrison-Zeldovich (scale- 
invariant, untilted) case is n = 1, which we adopt for our standard case. 

In contrast to the "standard" case in SHF, here we adopt the more con- 
servative choice (larger error bars) for the observed value of F, 

Fo = 0.25 ± 0.05 (8) 

(cf. Peacock and Dodds, 1994; Maddox, Efstathiou, and Sutherland, 1996). 
This is based on the galaxy correlation function, and it assumes that light 
traces mass. Equations (0) and (P) imply, very roughly, that Vly[h ~ 0.25. 

The shape-parameter constraint is in a sense the least robust of the con- 
straints we have discussed since it is not part of the basic Friedmann model. 
Rather, it depends on a theory for the primordial fluctuations and how they 
evolve. If the Friedmann cosmology were threatened by this constraint, we 
believe that those who model large-scale structure would find a way to dis- 
card it. Therefore SHF have also explored the consequences of removing 
this constraint and replacing it with one on derived from the M/L ratio 
in clusters and the luminosity density of the universe (Carlberg, Yee, and 
Ellingson, 1997) or one based on studies of large-scale flows around voids 
(Dekel and Rees, 1994; cf. Dekel, 1997). 



F VLyih exp —Ob — (/i/0.5) 



^/2(J7b/J1m)1 -0.32 (n-^ - 1) 



(7) 
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Figure 1. 68% (shaded) and 95% (dotted) confidence regions ("CRs") in tlie {Ho,Q.m) 
plane for CDM models with our four standard constraints. The CRs are closed curves. 
Individual constraints in this plane are also shown schematically. 

4. CDM Models: Results 
4.1. Standard Constraints 

For our standard case we have four observational constraints: ho = 0.70 it 
0.15, to = 14t2 Gyr, = fch^^'^ = 0.060 ±0.006, and To = 0.25 ±0.05. For 
this standard case we assume n = 1 and T = 0.9. The results for our three 
cosmological parameters are displayed in Figures 1 and 2 where the 68% and 
95% confidence regions ("CRs") are shown. Also shown are the projected 
CRs obtained by computing for single observables alone, or for pairs of 
observables. These are not true CRs but are intended to guide the reader 
in understanding how the various constraints influence the closed contours 
which show our quantitative results. 
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Figure 2. Same as Figure 1, but in the (Hoj'qio) plane. Individual and paired constraints 
are also shown schematically. 



Our best-fit values for the three key cosmological parameters are: r/io = 
8.2tli, Om = OAStlil and h = 0.58 ±0.22. Although the condition Quh « 
0.25 poses some threat to the SCDM (J^m = 1) model, Figure 1 shows 
that this threat is far from acute given our more accurate form of the T 
constraint in equation (0), as long as the error on ho is large (0.15) and BBN 
constraints are discarded. The exponential term in equation (0) becomes 
significant because the /g constraint forces f^B to increase with allowing 
the product ^luh to exceed 0.25. This has been noted before (White et 
al, 1996; Lineweaver et al, 1997). The SCDM model with Qm = 1 and 
h ~ 0.45 is acceptable but the high value for the baryon-to-photon ratio, 
7710 ~ 13, is in conflict with the inferred primordial abundances of all the 
light nuclides. Note that, although the uncertainties are large, low values of 
7/10 are disfavored (see Figure 2). 
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If we add the Dekel-Rees estimate of J^m (^o ^ 0.4), the five-constraint 
fit favors somewhat higher values of f^M and r/io and shghtly lower values 
of h. In contrast, if instead we include the cluster estimate (Oq = 0.2 it 0.1; 
Carlberg, 1997; cf. Carlberg, Yee, and Ellingson, 1997), we find a barely 
acceptable fit (Xmin ~ ^-0 ^oi 2 DOF, 92% CL), which favors lower values 
of and r^io and slightly higher values of h. 



4.2. Variations 



Tilt in the primordial spectrum has been investigated in many papers (Liddle 
et al, 1996a,b; White et al, 1996; Kolatt and Dekel, 1997; White and Silk, 
1996; Liddle and Viana, 1996). We considered the effect of a moderate "red 
tilt" {n = 0.8 instead of n = 1). This has the effect (see equation |^ of 
raising slightly the 68% and 95% contours in Figure 1. With this tilt the 
r constraint favors higher Om, so that the SCDM model is allowed for h 
up to nearly 0.5. The favored likelihood range for rjiQ is now also higher, 
though T^io ~ 7 is still allowed. However, the higher allowed range for r] does 
threaten the consistency of BBN. Conversely, a "blue" tilt, n > 1 (Hancock 
et al., 1994), would move the CR downward and allow models with J7m < 0.3 
at high h. 

Changing to a gas enhancement factor T = 1.3 (modest positive enhance- 
ment of gas in clusters) instead of 0.9 does not change the contours in Figure 
1 by much since F is only weakly coupled to VL-q through the exponential 
term in F. Although the effect is to lower the contours in Figure 1 slightly 
and to move downward the acceptable range for t/iq, t/io < 4 is still exclud- 
ed, disfavoring the low D abundance inferred from some QSO absorbers and 
favoring a higher helium abundance than is revealed by the Hii-region data. 

The possibility that the fraction of cluster mass in baryons in galaxies, 
isolated stars, and machos (/gal) might be larger - even much larger - than 
is implied by equation (|5|) would affect the CRs in much the same way as a 
small T, favoring even higher values of f^M and r/ig. 

The F constraint is crucial for our standard results favoring high and 
high r}. If, for example, we drop the F constraint and in its place use the 
cluster estimate Vto = 0.2 it 0.1, low Vty[ and low r/ are now favored (see 
SHF). 

The acceptability - or not - of the SCDM model depends crucially on 
the choice of Hubble parameter. SHF have experimented with replacing the 
standard constraint on Hq with = 0.70 ± 0.07. Now, the SCDM model is 
strongly excluded. 
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4.3. ACDM MODELS 



SHF have also considered models with nonzero A, limiting their investigation 
to the popular flat {k = 0) "ACDM" models with Q,\ = 1 — Q,m, where = 
A/ (3i?o )• For these models there are still only three free parameters and the 
four constraints discussed earlier are still in force, except that the product 
of the age and the Hubble parameter is a different function of 0,^ = 1 ~ ^A- 
t = 9.78 f {flu', k = 0) Gyr [Carroll, Press, and Turner, 1992, equation 
(17)]. For a given < 1, the age is longer for the flat {k = 0) model than 
for the A = model. The results differ very little from those in Figures 1 
and 2. The longer ages do allow the CRs to slide farther down toward large 
h and small J7m- Because of the longer ages at low Um (high ^^a), ^^o from 
clusters can now be accepted as a fifth constraint. In this case (see SHF) 
large Q,m and small h are now excluded while r/io > 4 is still favored strongly. 



5. Conclusions 

If BBN constraints on the baryon density are removed (or relaxed), the 
interaction among the shape-parameter (F) constraint, the cluster baryon 
fraction {fc) constraint, and the value of r/io assumes critical importance. 
These constraints still permit a flat CDM model, but only as long as h < 
0.5 is allowed by observations of h. The /g constraint means that large 
Qyi implies fairly large Ob- Therefore the exponential term in F becomes 
important allowing I^m = 1 to satisfy the F constraint. However, values of 
Vio ~ 8 — 15 are required (see Figures 1 and 2). The best-fit SCDM model 
has h !=a 0.45 and rjio 13, which is grossly inconsistent with the predictions 
of BBN and the observed abundances of D, ^Hc, and '^Li. For /;, > 0.5 a fit 
to SCDM is no longer possible. The SCDM model is severely challenged. 

The F and age constraints also challenge low-density CDM models. The 
F constraint permits flM < 0.4 only for high h, while the age constraint 
forbids high h, so Qm ^ 0.4 is required. The bound J7m ^ 0.4 conflicts 
with the added cluster constraint = 0.2 ±0.1 at the 98% CL, suggesting 
strongly that there is additional mass not traced by light. 

Although a few plausible variations on the CDM models do not affect the 
constraints very much, removing the F constraint would have a dramatic 
effect. Both high and low values of Om would then be permitted. The F 
constraint plays a crucial role in our analysis. 

At either low or high density, the situation remains about the same for 
the ACDM models. Because the ages are longer, we can tolerate « 0.3 
for h = 0.85. The ACDM model therefore accepts more easily the added 
constraint = 0.2 it 0.1. Improved future constraints on Q\ will come into 
play here. 
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Having bounded the baryon density using data independent of constraints 
from BBN, we may explore the consequences for the hght-element abun- 
dances. In general, our fits favor large values of r^io 6) over small values 
2). While such large values of the baryon density are consistent with 
estimates from the Ly-a forest, they do create some tension for BBN. For 
deuterium there is no problem, since for rjiQ ^ 6 the BBN-predicted abun- 
dance, (D/H)p ^ 3 X 10~^ (^c), is entirely consistent with the low abun- 
dance inferred for some of the observed QSO absorbers (Tytler et al, 1996; 
Buries and Tytler, 1996). Similarly, the BBN-predicted lithium abundance, 
(Li/H)p ^ 2.5 X 10^^'^ (^c), is consistent with the observed surface lithium 
abundances in the old, metal-poor stars (allowing, perhaps, some minimal 
destruction or dilution of the prestellar lithium). However, the real challenge 
comes from ^He where the BBN prediction for tjiq ^ 6, Yp ^ 0.248 {2a), is 
to be contrasted with the Hii-region data which suggest Yp ^ 0.238 (OS, 
OSS). 
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